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Abstract. We study the bound-state solutions of vanishing angular momentum in a 
quaternionic spherical square-well potential of finite depth. As in the standard quan- 
tum mechanics, such solutions occur for discrete values of energies. At first glance, it 
seems that the continuity conditions impose a very restrictive constraint on the energy 
eigenvalues and, consequently, no bound states was expected for energy values below 
the pure quaternionic potential. Nevertheless, a careful analysis shows that pure quater- 
nionic potentials do not remove bound states. It is also interesting to compare these 
new solutions with the bound state solutions of the trial-complex potential. The study 
presented in this paper represents a preliminary step towards a full understanding of the 
role that quaternionic potentials could play in quantum mechanics. Of particular interest 
for the authors it is the analysis of confined wave packets and tunnelling times in this 
new formulation of quantum theory. 

PACS. 02.30. Tb - 03.65. Ca 

I. INTRODUCTION 



In this introductory section, wc briefly recall the basic kinematic and dynamical framework of quater- 
nionic quantum mechanics. Obviously, it is not our purpose to give a complete and satisfactory 
introduction to the subject. For this and for further details on the mathematical formalism used in 
approaching a systematic development of quaternionic quantum mechanics, covering an axiomatic 
mathematical introduction to quaternionic quantum field and an interesting discussion on the role 
that quaternions could play in searching for new physics, we refer the reader to the excellent book of 
Adlcr [1]. 

In last years, the study of quaternionic quantum mechanics shed new light on the mathematical 
structures used in quantum mechanics. In particular, a better understanding of the eigenvalue prob- 
lem for real/complex and quaternionic linear differential operators [2,3] allowed the solution of the 
Schrodinger equation in presence of quaternionic potentials [4] without (unnecessary) complex trans- 
lations [5,6]. Nonetheless, quaternionic quantum mechanics is still in the early stage of development 
and what is still lacking is an interpretation of quaternionic potentials and proposals of experimental 
tests for quaternionic deviations of standard quantum mechanics. According to the nonrelativistic 
quaternionic scattering theory (developed in detail by Adler [1]) and the recent work on quater- 
nionic tunneling effect [7], experiments to detect left/right transmission asymmetry in presence of 
quaternionic barriers should involve CP-violating (asymmetric) potentials. Consequently, the earliest 
proposals for tests of quaternionic quantum mechanics based on neutron-optical experiments gave 
(as expected in virtue of the study on quaternionic scattering theory) a null result [8-10]. It is thus 
important to look for other possible signatures of violation of the standard quantum mechanics due 
to quaternionic potentials. The main scope of this paper is to solve the quaternionic Schrodinger 
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equation in presence of binding potentials. This could represent a potential (and not yet explored) 
research field to determine if and where quaternionic potentials could be seen. 



II. MATHEMATICAL TOOLS IN QUATERNIONIC QUANTUM MECHANICS 

To make our exposition self-contained and to facilitate access to the individual topics, the sections 
are rendered self-contained and, when necessary, we repeat the relevant material recently appeared in 
litterature. Before to begin the study of quaternionic binding potentials, we summarize our notation 
and introduce some basic concepts of quaternionic quantum mechanics. 

• Right multiplication by quaternionic scalars. 

In close analogy with the familiar complex quantum mechanics, the states of quaternionic quantum 
mechanics will be described by wave functions belong to an abstract vector space. Due to the non- 
commutative nature of the quaternionic multiplication, we have to specify whether the quaternionic 
Hilbert space, Vh, is linear under right or left multiplication by quaternionic scalars. Following the 
Adler book [1], we adopt the convention of linearity from the right for the quaternionic Hilbert space. 
Thus, if $„_i(r, and ^i^^if-it) belong to Vh and qi and ga arc quaternionic scalars, then 

$H,i(r,i)gi -l-$„,2(r,t)g2 e Vh . (1) 

• Anti-self adjoint operators and Schrodinger equation. 

The time evolution of <i>e(r,t) is governed by the quaternionic Scrodinger equation: 

9t$H(r-,t) = A(r,0$„(r,t) , (2) 

where 

A(r-,t) = [^|^V^-^K,a(^■,^)-JK,2(r,^)-fcT4,3(r,^)] /h 

is the anti-self-adjoint operator associated with the total energy of the system. The introduction 
of a pure quaternionic potential, jBc(r,t), does not modify the transition probability. There is an 
important difference between complex and quaternionic quantum mechanics. In complex quantum 
mechanics, we can trivially relate any anti-self-adjoint operator Ac(r, t) to the Hamiltonian observable 
TLcir^t) by the factor i, namely 

nc{r,t) =ihAc{r,t) . 

This is obviously not true in quaternionic quantum mechanics. In fact, the left-action of the complex 
imaginary unit i on j Bcir,t) still gives a pure quaternionic potential, kAc{r,t). 

• Stationary state and eigenvalue problem 

In the presence of time-independent potentials, the quaternionic Schrodinger equation (0) becomes 

at$H(r-,t) = A,(r)$H(r-,t) 

*i&^'-*"^«,iW"j''^«.2(r)-fcl4,3(r-)] Mr,t)/h (3) 

The method of separation of variables is still applicable if we set 

$n(r,i) = fMMxmit) ■ 

Due to the linearity under right multiplication by quaternionic scalars of Awi{r), the right position 
of the time-dependent quaternionic function Xu{t) plays a fundamental role in separating the initial 
differential equation in two equations containing respectively the time and space variables. Explicitly, 
we find 

^XhWXh^W^V'h^W f*i^V^-zK,iW-jK,2(r)-fcK,3(r)l Mr) ■ 
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This equality is only possible if each of these functions is a constant, which we shall set equal to A e H. 
Thus, the right linearity of the quaternionic Hilbert space Vh implies a right eigenvalue problem for 
the anti-self adjoint operator Amir), 

A,-,{r) ip,,{r) = ip,,{r) X/h . (4) 

The need to use right eigenvalue in physical problems can also be demonstrate by simple mathematical 
arguments [2]. 

• Number systems and amplitudes of probabilities. 

To give a satisfactory probability interpretation, amplitudes of probability must be defined in asso- 
ciative division algebras. Amplitudes of probabilities defined in non-division algebras fails to satisfy 
the requirement that in the absence of quantum interference effects, probability amplitude superposi- 
tion should reduce to probability superposition. The associative law of multiplication (which fails for 
the octonions) is needed to satisfy the completeness formula and to guarantee that the Schodinger 
anti-self-adjoint operator leaves invariant the inner product. The important point to note here is that 
we can still use vectors in Clifford algebraic or octonionic Hilbert space. Nevertheless, amplitudes 
of probability must be defined in C or H. For example, we can formulate a consistent complexified 
quaternionic [11] or octonionic [12] quantum mechanics by adopting complex inner products [13], i.e. 
complex projections of complexified quaternionic or octonionic inner product. The use of complex 
inner product represents a fundamental tool in applying a Clifford algebra formalism to physics and 
plays a fundamental role in looking for geometric interpretation of the algebraic structures in rela- 
tivistic equations [14-17] and gauge theories [18, 19]. The choice of quaternionic inner product seems 
to be the best adapted to investigate deviations from the standard complex theory. In this paper, 
we shall work with quaternionic inner products. Thus, with each pair of elements of the quaternionic 
Hilbert space Vh, we associate a quaternionic inner product defined by 

(</Jh,'0h> := J dr lp^^{r) ip^r) e H, (5) 

where ^j,(r) denotes the quaternionic conjugate of 

(fnir) = (pwL,o{r) +iiptt,iir) + j ipM,2{r) + k ipt,,3{r) , 

i.e. 



• Energy eigenvalues. 

The anti-hermiticity of Aa{r), 

{ ipa , Amlpa) = - { Aa ifia , V'h ) 

implies A = — A. Consequently, observing that A is related to the total energy of the system and setting 

A = - {iE, + jE, + kE,) , 

the quaternionic energy eigenvalue equation Q becomes 

Ai(t-) Mr) = - Mr) {iE, + jE^ + kE,) /h. (6) 

In the complex limit, the previous equation reduces to 

Hc{r) ifc{r) = ihAc{r,t) (pc{r) = E^ ifc{r) 

Eigenvectors of anti-self-adjoint quaternionic operator divide into mutually orthogonal eigenclasses. 
Each eigenclass corresponds to a ray of physically equivalent states. Eigenvectors within each eigenclass 
are not orthogonal and have eigenvalues related by the automorphism transformation 

u {i Xi + j X2 + k A3) u , 
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with u unitary quaternion. By choosing a particular automorphism transformation, i.e. 

A9 ^Aq 



2j\l+\l+\l 



1+] 



we set a particular representative ray within each eigenclass whose eigenvalue is complex, 



A = i-y/AJ +A3+A3 . 
Hence, within each eigenclass of Auir) eigenvectors, there is a ray 

for which 

A,(r)*H(r) = -■^^{r)iE /h . 
Thus, for stationary states, the quaternionic Schrodinger equation can be rewritten as follows 

|i - j K,,(r) - j 14,3 (r) - k K,3(r) ] *„(r) = - 4'„(r) i E . 



(7) 
(8) 



III. SPHERICALLY SYMMETRIC POTENTIALS 

The central force problem in which V^^i{r) = occupies a prominent place in (complex) quantum 

mechanics because many important applications involve potentials that are at least approximately 
spherically symmetric. In this section, we extend this symmetry to the quaternionic part of the 
potential, 

Vj,^2{r) = 14,2 and K.sl'^) = K,3('') • 
Since the quaternionic potential 

V,,{r)=i (r) +jVu4r)+kVu,, (r) , 

depends only on the distance r of the particle from the origin, spherical coordinates are the best 
adapted to the problem. Thus, we express the Laplace operator in Ati{r) in spherical coordinates 
(r, 0, Lp) by the well-known formula 
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Since the operator depends only on 9 and tp, by using the (complex) spherical harmonic functions 

*H(r-)=i?„(r)y,""(0,^) , 
we can reduce the Schrodinger equation (jSj) to the following quaternionic radial equation 



2m 



dr^ r dr 



where 



(9) 



2m 



By the change in the function 



i?H(r-) =C/H(r)/r , 
we obtain for Uti{r) the following differential equation 

Zm dr'' 



U,,{r)^~U^{r)iE 



(10) 
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lll-A. SPHERICAL POTENTIAL TRAP AND QUATERNIONIC EIGENFUNCTIONS 

To find the energy levels of a particle with angular momentum / = in presence of a quaternionic 
spherical potential trap 



r < a (region I) , 

iVi+jV2 + kV3 r>a (region II) , 

we have substitute such a potential in Ea. (|10|l . This yields the two ordinary differential equations 

zC^i,H(r) = -C/:,„(r)z2m^/ft% (11) 

and 

iU^Ar)-['2m{iV^+jV, + kV,)/ff]UuAr) ^ - UnAr) ^ ^rnE / , (12) 

The solution in region is easily obtained by observing that to the standard complex sine/cosine solution 
we have to add a pure quaternionic hyperbolic sine/cosine solution, 

Ui.uir) = sin(er) Ui +cos(er) f3i + j sinh(er) 7i + j cosh(er) 5i , 



where e = ^j2mE /h? and ai , /?i , 71 , 5i are constant complex coefficients to be determined by bound- 
ary conditions and continuity constraints. To guarantee the correct behavior of the solution at the 
origin, the wave-function must satisfy f7i ,,(0) = 0, thus we have to impose the following constraint 
/3i = (5i =0. Consequently, the solution inside the well (r < a) is given by 

UiA'r) = sin(er) Qi + j sinh(er) 71 . (13) 

The discussion in region II is a little more complicated (see the appendix A for the explicit calculations) 
and provides the following solution 

UiiA''') = {1+ jw) { sinh[z^_ r] + cosh[v_ r] 13^} + {z + j) { sinh[zy+ r] 7^ + cosh[i^+ ''l '^2 } , 

where 



^ (v,±JE''- - 1 , w^-i ^^=^ and 



V2+iVi 



e+^e'-v'^-vI e + ^e'-v'^-v. 



Observing that, independently of the value of i?, the real part of i/+ is never null, the condition that 
the wave function does not diverge at r = -l-oo requires a first constraint i.e. 72 = —b^ (without loss 
of generality we can choose Re[;^-t] > 0). Thus, the solution in region II can be rewritten as 

t^ii,n('') = (1 + J w) [sinh(z^_ r) a.^ + cosh(i'_ r) p.^] + {z + j) cxp[-i^+ r] 62 , 



If the total energy E is greater than \/vf + V^+V^ , we have the free particle case. In fact, for 
E > \J Vl + V2 + , iy_ is a complex imaginary number and no constraint exists for the coefficient 
a2 and /Sj • The continuity relations in r = a for the wave- functions and their gradients can be satisfied 
for all values of the energy and the energy levels are therefore not quantized. 

The bound particle case corresponds to E < \/V^ + + . In this case, the condition that the 
wave function docs not diverge at r = +00 gives an additional constraint on the coefhcients aa and 
(32, i.e. = — /Ja- Thus, we have a spatially damped solution in region II given by 

Uu,u (r; E < ^JVt + V'i + V'i ) = (1 + j w) exp[-z/_ r] /S^ + (z + j) exp[-j/+ r] ^2 ■ (14) 

It remains to adjust the coefficients «!, 71, /?2 and so that the amplitudes and gradients at r = a 
are continuous. This will imply a quantization for the energy levels. 
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lll-B. QUATERNIONIC BOUND STATES 

The continuity conditions 

C/i,H(a) = J7n,n(a) and [/,,H(a) = J7ii,n(a) 

imply 

sin(ea)ai + j sinh(ea)7i — cxp[— i^_a] /Ja + z cxp[— i^^a] 62+ 

j [w exp[-i^_a] /Ja + exp[-j/+a] 82] , 

— e [cos{ea) +j cosh(ea)7i] = exp[— i^_a] (32 + exp[— i^+a] ^2+ 

j [wiy^ cxp[— j/_a] ^2 + 1^+ cxp[— i^+a] ^2 ] ■ 

Separating the complex from the pure quaternionic part in the previous equations and eliminating 
the coefficients and 71, we obtain the following matrix equation 



i'_ tan(ea) + e z [i/^ tan(ea) + e]\ / exp[— a] (32 
w [v_ tanh(ea) + e] tanh(ea) + e j I cxp[— a] ^2 

This matrix equation has a non-trivial solution only if the determinant of 

tan(e a) + e z [u^ tan(e a) + e ] 
w [ v_ tanh(e a) + e ] Vj^ tanh(e a) + e 

vanishes, i.e. 



z w 



[ v_ tan(e a) + e ] [ tanh(e a) + e ] 
[ v_ tanh(e a) + e ] [v^ tan(e a) + e ] 
This equation can be rewritten in a more convenient form, i.e. 

[u^ ~ zw v_) tanh(e a) + ( 1 — 2 w ) e 



tan(e a) ~ 







(15) 



(16) 



{ \ — zw) tanh(e a) + { \ — z w f+jv^ ) e 



tanh(e a) 



1 - 



1 - 



tanh(e a) 



(,2 + y^T^) 



(17) 



where 



and 



Let us first examine the case: \/V^~^-V^ < E < \JVi + V2 + yl ■ Observing that this condition on 
the energy eigenvalues implies Vq < t and > (e^* — i^qY^'^, it can be immediately seen that Ea. (|17|l 
represents a real equation which generalizes (to the "quaternionic" case) the well-known equation 
obtained in "complex" quantum mechanics, i.e. 



tan(ea) = /(e; i/c, 0) 



(18) 



One may ask wheter this is still true for energy values below y/vT+V^. For E < ^V^ Vl , it can be 
immediately seen that v± and zw are complex and this seems to be too restrictive. Consequently, no 

bound states should be expected for E < \/Vl + . Nevertheless, by simple algebraic manipulations 
(see appendix B for the explicit calculations), it can be shown that Ea. ll7() still represents a real 
equation and bound states could exist below the pure quaternionic potential. Graphical solutions of 
equation (|17|l give the energies of the bound states of a particle in a quaternionic spherically symmetric 
potential trap, see Fig.l and Fig. 2. It is also instructive to compare the quaternionic bound states 
with the bound states of the trial-complex potential y/V^ +V2 +V^. 
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IV CONCLUSIONS 

This work was intended as an attempt at motivating the study of quaternionic formulations of quantum 
mechanics and at looking for deviations from the standard complex theory. The analysis presented 
in this paper represents a preliminary discussion on quaternionic bound states and touched only a 
few aspects of the theory. In particular, no attempt has been made here to develop a quaternionic 
perturbation theory or to study the behavior of confined quaternionic wave packets. These topics 
exceed the scope of this paper but surely represent important points to be investigated in view of a 
complete understanding of how and where quaternionic potentials could be seen and if they really 
exist in nature. 
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Appendix A 

To calculate the radial solution in region II, we have to solve a second order differential equation with 
quaternionic constant coefficients, 

i ^ tV, ~ - kvA UnAr) = - UuAr) i E . (19) 

2m ar^ J 

We refer the reader to [4, 20] for a detailed analysis of quaternionic differential operators. In this 
appendix, we will touch only a few aspect of the theory, by restricting our attention to differential 
operators with quaternionic constant coefficients. A (right-complex linear) solution of Ea. H19(l can be 
written in terms of (left-acting) quaternionic (g) and complex [v) coefhcients, i.e. 

q cxp[j/ r] . 

To determine this coefficients, let us apply to Ea. H19|) the anti-Hermitian operator 

d' 

Zm ar-' 

In this way we obtain 

y d" _ _ d 



2m / dr'i 



2 V^i ^ - {V," + V;' + ^3') I UnAr) =-E" UnAr) , (20) 

which represents a real differential equation. Consequently, the quaternionic factor q can be factorized 
and the complex coefficient v is calculated once solved the following algebraic equation 



— „^-2 — V,iy' + V,+V,+V,-E^0. (21) 
2m / 2m 



The solutions for the complex coefficients v are given by 



1/1,2 = ±v- and = ± i/+ , where = y'^ (v, ± A - V,^ - ) . (22) 

Coming back to Eq.linil, for E ^ y/WTvl and E ^ V^i' + V^' + we then find four (right- 
complex) linear independent solutions [4], i.e. 

(j_ exp[— r] , q_ exp[i/_ r] , q^ exp[— i/+r] and q^ exp[i/+r] . 

To find the quaternionic factor, q_ (associated to exp[± v^r]) we set q_ = \ + j w {w Cz C). This choice 
is possible due to the right-complex linearity of the quaternionic differential equation (|19|) . It follows 
immediately that 



^{l+jw) (v,-^/E^^v[~^^ -itV,+jV, + kV,){l+jw)^~{l+jw)'. 
previous equation implies 
(^E - ^Je^' -V,' -V^^ + (X4 + i V,) «; = , 



\iE 

It is easy to check that the previous equation implies 

ilE 



i [E + ^^E~ - < j w - (y^ - i V;) = . (23) 

Finally, 

. E - AE' -V^ -V: V.-iV^ .... 

w = — I = —i , (24) 

V^+iV, E + Ae' - vl - vl 

To find the quaternionic factor, q^ (associated to exp[±i^+r]) we choose (7+ = z -t- j (z G C). A 
calculation similar to previous one gives 



, ^ E-AE'-y:-y: _, ^ + »^ 

v,-iv, E + Ae' -v:^~ ■ ^ ' 
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Appendix B 



In this appendix, we are interested to prove that Ea. (|17|l for E < y^V^ + still represents a real 
equation. Before to proceed with the proof, let us observe that 



zw = exp 



-2 i arctan 



E 



and 



v± 



^ + ^v: + v: + v:-E-'±tj^ ^^v: +v: +v: -e^ -v. 

It is natural to rewrite Ea.()17|l as follows 

{v^ — zw tanh(e a) + ( 1 — z w ) e 



tanle a) 



v_v+ [1 — zw) tanh(e a) + {v_ — zw v+) e 
NuM NuM Den 



Den 



|Den|- 



(26) 



The important point to note here is that 

= ~ and \zw\ = 1 . 

We can now proceed with the proof, 



NuM Den = [{v^ — zwiy^) tanh(e a) + {1 - zw) e][\v_\'^ (1 — zUJ) tanh(ea) + ( i/^ —'zwv_ ) e] 
= {v_ — zw v_ + ~ — zw~ ) 1^ tanh^(e a) + 

[ + ~^ — {zw -\- zUJ) li'- 1^ + (1 — 2 w) (1 — Tw) I'y'-l^ ] e tanh(e a) + 

( — z w v_ + ~ — zw vZ ) 
e M . (27) 
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Fig. 1 



tan(e a) 

f{ea- 5 7r, 0) 

f{€.ci\ 5 TT, 2.5 tt) 

/ (ea; 1.0625^/^ - 5 tt, 0) 



bound states [e a] 



-3 - 



o 

2.9525 
5.8983 
8.8280 
11.7238 
14.5273 



2.9561 
5.9048 
8.8374 
11.7378 
14.5541 



O 

2.9552 
5.9040 
8.8376 
11.7391 
14.5578 



f 



o 



8 



12 



e a 



Fig. 1. Graphical solutions of equation 11711 . The circle, bullet and diamond symbols respectively represent 
the vanishing angular momentum bound state solutions of the radial Schrodinger equation for a complex 

( i Vi ), quaternionic ( i Vi + j V2 + V3 ), and trial-complex {yi\/vt + + j square-well potential. For 
the convenience of the reader we also list the n umerical values of the bound energies. The plots refer to the 
following potentials: a^/2niVi/W = 5 tt and a \J 2m \/V^ -f / h? — 2.5 vr. 
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-3 - 




O 



8 



12 



16 



e a 



Fig. 2. Graphical solutions of equation 11711 . The circle, bullet and diamond symbols respectively represent 
the vanishing angular momentum bound state solutions of the radial Schrodinger equation for a complex 

(iVi), quaternionic {iVi + j V2 + kVa), and trial-complex ^ i y/ + V2 + ^ square-well potential. For 
the convenience of the reader we also list the n umerical values of the bound energies. The plots refer to the 
following potentials: a\j2m Vi/h? = Stt and a \J 2m y/V^ + Vi / ft^ = Svr. 



